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Thermal boundary conductance at a metal-dieletric interface is a quantity of prime importance 
for heat management at the nanoscale. While the boundary conductance is usually ascribed to the 
coupling between metal phonons and dielectric phonons, in this work we examine the influence of 
a direct coupling between the metal electrons and the dielectric phonons. The effect of electron- 
phonon processes is generally believed to be resistive, and tends to decrease the overall thermal 
boundary conductance as compared to the phonon-phonon conductance a p . Here, we find that the 
effect of a direct coupling a e is to enhance the effective thermal conductance, between the metal 
and the dielectric. Resistive effects turn out to be important only for thin films of metals having 
a low electron-phonon coupling strength. Two approaches are explored to reach these conclusions. 

First, we present an analytical solution of the two-temperature model to compute the effective 
conductance which account for all the relevant energy channels, as a function of a e , er p and the 
electron-phonon coupling factor G. Second, we use numerical resolution to examine the influence of 
<7 e on two realistic cases: gold film on silicon or silica substrates. We point out the implications for 
the interpretation of time-resolved thermoreflectance experiments. 


I. INTRODUCTION 

Interfacial thermal transport is of major importance in heat application management. Indeed, ohmic contacts are 
ubiquitous in microelectronics and experience increasing levels of loads. As a result, a fundamental understanding 
of the cooling kinetics of heated metals or semiconductor at ohmic contacts is a prerequisite to design electronic 
components that can support large heat flux densities [jJ. Thermal relaxation of microelectronics components is 
primarily governed by the boundary conductance, since as the density of interfaces increases, thermal losses in the 
bulk become negligible compared to interfacial resistance. While the boundary resistance involves a priori the energy 
transfer between metal electrons or metal phonons, and the phonons in the substrate, most of the models developed 
so far consider only the latter channel of energy. These include the acoustic mismatch model (AMM), where phonon 
transmission at an interface is assumed to be controlled by the difference in the acoustic impedances of the two 
media [2], and the diffuse mismatch model (DMM) which posits that interfacial roughness destroys any correlation 
between incident and reflected phonons [|j. 

The boundary resistance of metal-diamond interfaces at room temperature has been first experimentally determined 
by Stoner and Maris in 1993 p|, using time-resolved thermoreflectance (TTR). This work has been followed by the 
characterization of various interfaces, extending the early measurements to higher temperatures M- More recently, 
TTR has been applied to investigate heat transfer across metal-graphite Q and metal-carbon nanotube interfaces 
as well 0- In parallel, theoretical progress has been achieved through molecular dynamics 00 which ignores 
electronic degrees of freedom, or treat them in the spirit of the two-temperature model [l3j]. A recurrent conclusion 
of this body of works is that the Kapitza conductances measured between two solids depart significantly from the 
classical AMM and DMM models, which often yield comparable values. For instance, the thermal conductance of 
one of the most studied system, the Pb-diamond interface, is found to be typically sevenfold higher than the DMM 
prediction. 

After the early room temperatures measurements of Stoner and Maris, several authors have proposed to relate this 
excess of conductance to electron-phonon scattering taking place across the interface EM3- Indeed, there are at 
least two mechanisms according to which electrons may heat up the substrate, as represented in Fig. [l] The first 
one has been extensively studied and involves an indirect coupling between the electrons and the substrate through 
phonon-phonon processes fl3l ITil I~ilil22[ . This mechanism is only operative when the electrons are not in equilibrium 
with the lattice. In this situation, the effect of the electron-phonon process is to introduce an additional resistance 
1/y/Gkp where G is the electron-phonon coupling factor, and k p the phononic thermal conductivity of the metal. 
Hence, it is found that electron-phonon processes deteriorate the interfacial energy transfer, as compared to the case 
where electronic degrees of freedom are neglected. 

There is another energy channel that has been recently unveiled from the experimental side (23l425| . and which 
has been previously predicted theoretically EKnunt This mechanism is a direct coupling of the electron-phonon 
coupling through the interface, characterized by an interfacial heat flux of the form cr e (T e — T s ), which depends on 
the temperature difference between the electrons and the phonons substrate, through the electron-phonon interfacial 
conductance a e . Throughout the decade 2000 and until recently, the latter electron-phonon transfer mechanism has 
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FIG. 1: (Color online) Schematic representation of the two energy transfer channels between metal electrons (e _ ) and substrate 
phonons (ph). The indirect coupling (blue) involves electron-phonon coupling in the bulk metal, characterized by G followed by 
transfer via the phonon-phonon conductance a p . This channel has been studied by Majumdar and Reddy [16}. The direct route 
studied in this work and represented in green, involves the electron-phonon conductance <r e , for which theoretical predictions 
exist [Til . FiTI . [It} • Diffusive transport represented by dots plays a role in both cases. 


been disregarded essentially because experiments where the metal is replaced by Bi a semi-metal and which conclude 
that the values of the interfacial conductance seems to be independent on the nature of the electronic transport 
properties of the film Q. There are, however, actually two hints that the direct electron-phonon channel could play a 
significant role. First, the above mentioned models for the conductance cr e predict values that are comparable or even 
largely exceed the typical values for the phonon-phonon conductance <r p . Second, recent transient thermoreflectance 
experiments on a Au film [ 23 ],[24], [26| concluded on the existence of an electron-substrate energy transfer, with estimate 
of a e in the range 100 to 1000 MW uU 2 K -1 , compatible with the theoretical predictions. The consequences of this 
direct channel of energy have not been discussed in the literature. 

In this article, we aim at a theoretical description of the competition between the three different channels represented 
in Fig. [I] We show that electron-phonon processes taking place in the bulk of the metal induce a small, commonly 
negligible resistance, while the effect of the direct coupling is to enhance interfacial heat transfer. Importantly, we 
find that the effect of a direct coupling persists over long time scales on the nanosecond range, even though the 
electrons and the phonons in the metal have been long allowed to equilibrate. Indeed, if the largest values of a e 
are experimentally reported in situations where electrons are out of equilibrium with the metal lattice [24[, still a 
non vanishing value a e > 100 MW m -2 K -1 is measured in case of negligible electron-phonon nonequilibrium in the 
metal [23j. This indicates that this channel of interfacial transfer should be operative in the TTR measurements of 
metal-dielectric interfaces, and may explain the large discrepancies observed between the DMM predictions and the 
TTR measurements. 

The article is organized as follows: In Sec. pi] we first motivate the study by calculating the interfacial electron- 
phonon conductance for different interfaces using Sergeev’s model for a e . Using analytical calculations on the two- 
temperature model, we examine in Sec. IIIII the influence of a finite <r e on the temperature decay of a heated metal. In 
Sec. IIV1 we illustrate the effect of a e for two specific systems: gold/silicon and gold/silica interfaces. We show that 
for electron-phonon conductance larger than the phonon-phonon conductance, it is actually the dominant mechanism 
governing the temperature decay. Our conclusions are given in Sec. [V] 


II. THEORETICAL ESTIMATES FOR ELECTRON-PHONON AND PHONON-PHONON 

CONDUCTANCES 

In this section, we motivate our study by estimating, for various metal/non-metal systems, the electron-phonon 
and phonon-phonon interfacial conductances. In both cases, we rely on theoretical predictions. 
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Electron-phonon conductance a e . There are three analytical models predicting; the electron-phonon interfacial 
conductance at a metal/diclectric interface. Huberman and Overhauser [14| first propose that metal electrons couple 
with joint vibrational modes at the interface leading to a significant enhancement of the energy transfer, as compared 
with direct phonon-phonon coupling. Later, Sergeev developed a field theory to predict the effect of the boundaries 
on the inelastic scattering of conduction electrons [l5[. Finally, Mahan derived a theory where metal electrons may 
transfer heat to the substrate through the image charges created by the vibrating ions of the non-metal 0 . This model 
should mainly apply to ionic crystals or polar semi-conductors, rather than to apolar non-metals. In both Huberman 
and Mahan’s models the conductance is predicted to be constant at high temperatures, while in Serg eev’s model, it 
is found to increase linearly with temperature [39| . in agreement with experimental observations |23| . Because it is 
the sole model to capture this temperature dependence, we chose to consider only Sergeev’s model in the following. 

As far as we know, there has been no attempt to quantify the values of the electron-phonon interfacial conductance 
predicted by Sergeev’s model for real cases at ambient temperatures. Here we determine the values of a e for different 
metal/dielectric interfaces. We do so by using input data from density function theory (DFT) calculations predicting 
the electron-phonon coupling of several noble metals. To proceed, we start with the expression 


cr e 


3 ?rh 7 u L + 2 {H±\ 
35((3)k B T e _ p \u T ) 


(1) 


where h is the reduced Planck constant, k B is Boltzmann constant, ( is Riemann zeta function, 7 is the Sommerfeld 
constant of the metal, u B and ut are respectively the metal longitudinal and transverse sound velocities and r e _ p is 
the electron-phonon energy relaxation time of the bulk metal. Because at high temperature 1/T e _ p ~ T e [28| while 
other factors have weak or no temperature dependence, one immediately obtains a e ~ T. Here, we estimate the 
electron-phonon energy relaxation time using its relation with the electron-phonon coupling factor Hill]: 

( 2 ) 

6 T e —p± e 

where m e is the effective mass of the electron, c s is the speed of sound and n e is the free electrons number density. 
Besides, one has: 


G _ 3h\{uj 2 ) 

lT =~C e - 7Tk B T e ’ 


( 3 ) 


where A is the electron-phonon coupling constant, so that the electron-phonon conductance Eq. [Qmay be cast as 


_ 54 7 2 ft 2 A(o; 2 ) 

° 33t: 2 k 2 B m e c s n e 

To calculate a e , we have used the values of A(w 2 ) from the DFT calculations of Lin et al. [28| for noble metals. The 
effective mass to* is taken equal to the electron mass m e , except for Platinum, for which we have used m* = 13 m e [30| . 
For Bismuth, we have calculated the value of A(w 2 ) from the thermodynamic analysis of Giret et al. [3j], H2], who 
reported an electron-phonon coupling factor G = 1.45 10 16 W m -2 K -1 . Table U gathers the values of the electron- 
phonon interfacial conductance for common metals and the semi-metal Bi. 
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metal 

n 

ul 

ut 

n e 

7 

A(c 2 > 

CT e 

Al 

100000 

6240 

3040 

18.1 

1.3 

185.9 

2.64 

Au 

97970 

3390 

1290 

5.9 

0.67 

23 + 4 3.08 

Bi 

46794 

1543 

1107 

14.1 

0.08 

229.3 

0.02 

Cr 

138269 

6980 

4100 

16.6 

1.57 

128 

6.2 

Pb 

55990 

2350 

970 

13.2 

2.93 

45 + 5 

17.8 

Pt 

110872 

4174 

1750 

3.2 

6.54 

142.5 

186 


TABLE I: Physical parameters for the metals: molar density n (molm -3 ), longitudinal and transverse sound velocities ul and 
ut (ms -1 ), free electron number density n e (10 28 m -3 ) , Sommerfeld’s constant 7 (mjmol -1 K -2 ), factor A(w 2 ) (meV 2 ), and 
electron-phonon conductance o e (100 MWm - 2 K -1 ) at room temperature. 
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Phonon-phonon conductance er p . Here we use the diffuse mismatch model Q, as we briefly discuss now. Quite 
generally, the phonon-phonon thermal boundary conductance is given by 




Tp - df( 

Vi iP hojDi iP (u>) 


eq 


(co,T) 


rir/2 


dT 


a p (oj, cos 6) cos 6d {cos 9)dw, 


(5) 


where the sum runs over the polarizations p , uf p is the cut-off frequency, V\, p is the mode group velocity, Di p (lo) is 
the density of states, / eq is the equilibrium Bose-Einstein distribution, and a p (cu,cosd) is the phonon transmission 
coefficient. If interfacial scattering is supposed to be diffuse, the transmission coefficient no longer depends on the 
angle of incidence, and from detailed balance it follows that 


E f 


a = 


E 


p v f 


E 


p vi 


( 6 ) 


where the density of states has been taken to be given by Debye’s model. Note that the transmission coefficient is 
assumed to be polarization independent [33]. To compute er p , we have used the acoustic properties of the different 
metals and dielectric provided in Tabs. [T] and HIl 


substrate 

n (molm d ) 

ul (ms Q 

ut (ms 1 ) 

A1 2 0 3 

38921 

10890 

6450 

AIN 

160345 

11120 

6267 

diamond 

292667 

17500 

12800 

GaN 

138269 

7781 

4427 

Si 

73214 

8970 

5332 

SiOa 

44167 

5950 

3740 


TABLE II: Physical parameters of the substrates: molar density n, longitudinal and transverse sound velocities ul and ut- 

Figure [2] gathers the electron-phonon and phonon-phonon conductances for the different systems considered. It is 
apparent that in most cases <r e is much larger than er p . This implies that for most of the metal/diclectric interfaces, 
tr e may not be neglected, and should be accounted for in the analysis of thermoreflectance data. Qualitatively, one 
expects that compared to the bare phonon-phonon conductance, the apparent thermal conductance of the systems 
discussed should be enhanced, as will be investigated quantitatively in the next sections. Besides, note that even 
Bismuth which is a semi-metal, and not considered a good metal is predicted to have a significant a e . This suggests 
that the relevant energy channel discussed here should be present not only in metals, but also in semi-metals despite 
the fact that the number of free electrons is relatively small. Hence, even materials which display poor electronic 
transport properties such as semi-metals, could be affected by electron-phonon interfacial coupling. This is because 
for these materials the electron-phonon relaxation time takes values slighty larger but comparable with those of good 
metals. 

It should be kept in mind however, that our estimates of the transport coefficient cr e comes from a model which 
may have some limitations. In particular, Sergeev’s prediction yields a conductance which depends only on the 
bulk properties of the metal, and not on the nature of the substrate. Also, it should be remarked that for gold 
systems, Sergeev’s model overestimates by a factor of two the electron-phonon conductance reported experimentally, 
as shown below in Fig. [3] Further work is clearly needed to develop models of interfacial electron-phonon coupling, 
which includes information from the interface. However, even if Sergeev’s model may overestimate cr e , a e remains 
comparable to cr p in many cases, suggesting that from those theoretical estimates, there is no reason to discard the 
electron-phonon conductance. Given the uncertainty of theoretical estimates for a e , and the lack of experimental 
measurements for most systems, we will not choose any specific value or model, but rather explore a wide range of 
values for a e . 


III. ANALYTICAL PREDICTIONS 

We derive in this section analytical expressions for the characteristic cooling time of a heated metal/dielectric 
interface, taking into account not only the phonon-phonon conductance, but also the electron-phonon conductance. 
The analytical solutions may allow to understand quantitatively the effect of a e on the effective interfacial conductance 
which account for the different channels, and analyze the limiting behaviours. 
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FIG. 2: (Color online) Phonon-phonon conductance <j p and electron-phonon conductance a e estimated from Sergeev [5j| and 
the DMM model respectively, at room temperature. The black solid line indicates the curve a e = a p . 



FIG. 3: (Color online) Temperature dependence of the electron-phonon interface conductance for gold/silica and gold/silicon, 
as reported experimentally [23|]. The solid lines are the values obtained from Eas. 1171 and 1181 The blue dashed line shows 
Sergeev’s predictions for gold, as given by eq. [4] 

A. Model and parameters 

The situation relevant to TTR experiments is schematically depicted in Fig. [31 a metal film with thickness h on a 
semi-infinite dielectric substrate [4(|. To keep calculations tractable, we assume that the problem is one-dimensional, 
with variation occurring only perpendicular to the interface (^-direction). In principle, the cooling kinetics should be 
described by the solution of the Boltzmann transport equation, but the long time regime on the order of 1000 ps that 
we discuss here is generally well captured by the two-temperature model where the electrons and phonons are each 
characterized by a temperature, subject to diffusion and relaxation [2l], |29j, [34J • The evolution equations are 

c e d t T e = k e d 2 xx T e - G(T e -Tp), 
c p d t T p = fc p d' xx T p + G(T e - T p ), 
c s dtT s = k s d xx T s . 
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FIG. 4: (Color online) Sketch of the system modeled: a thin metallic film on a dielectric substrate, and irradiated by an ultrafast 
laser pulse. T e , T p are the electron and phonon temperatures in the metal, T s is the phonon temperature in the substrate. 


Here, we use subscripts i = e, p, and s for the metal electrons, the metal phonons and the substrate phonons 
respectively, Ci and k.-, are the heat capacities and thermal conductivities, G is the coupling factor. Unless otherwise 
mentioned, we will consider that all coefficients are constant, and in particular independent of temperature. All 
quantities are given in SI units, unless indicated otherwise. As regards the boundary conditions, we assume a zero 
energy flux at the air-metal interface, while the metal-dieletric interface is characterized by two boundary conductances, 

for x = 0, —k e d x T e = 
for x = h, —k s d x T s = 


Here a p and a e are the phonon-phonon and electron-phonon conductances respectively and the second equation 
expresses the continuity of the energy flux across the interface. Finally, the initial conditions are 

for f = 0, T e = T\x), T p = T s = 0. (7) 

T e , T p and T s are thus the differences with respect to the temperature T 0 of the system prior to any laser heating [4l|. 
The initial temperature profile T l is exponential with decay length 5, the penetration depth of the laser in the metal 
film: 


-fc p d x T p = 0, 

cr e (T e - T s ) + cr p (T p - T s ), 


T-(x) = T*(0) e-/‘, T{ 0) = S{1 _ h e . kll y 

The normalization constant T®(0) is such that the average temperature increase in the film is unity [42.]. Below, we 
focus on the electron temperature at the surface of the film, and its long-time behavior [43[ . 

Tabic IIIII lists the various quantities needed in the problem, from which one can construct the relaxation, diffusion 
and conductance times. The relaxation time Ti = Ci/G is related to the electron-phonon coupling in the bulk of 
the metal film, = h 2 /on is the typical time for diffusion over the film, and r C j = Cj/i/cq characterizes the cooling 
of the carriers induced by the interfacial conductances. Taking ratios between those characteristic times leads to 
dimensionless parameters: the R number compares the electron and phonon relaxation times, the Biot numbers 
Bi =0 ,p compare the diffusion and conductance times and the M, =e , p numbers compare the relaxation and diffusion 
times. Below we focus on the influence of the two interface conductances, with <r e taken in the broad range 10 5 — 10 11 , 
and dp in the range 10 7 — 10 9 W m -2 K -1 , which should encompass most cases of experimental relevance. All other 
parameters are kept fixed to a value representative of real systems, those default values are also given in Tab. IIIII 


B. Influence of the electron-phonon conductance o e 

We now investigate how, within the two-temperature model, a finite electron-phonon conductance a e affects the 
temperature decay. We first present a simplified model, whose analytical solution is straightforward. 
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Quantity 

Symbol 

Definition 

Value 

Unit 

film thickness 

h 

- 

100 

nm 

penetration depth 

S 

- 

20 

nm 

heat capacities 

Ce 5 Cp, Cs 

- 

10 4 , 10®, 10® 

J m -3 K _1 

thermal conductivities k e , k p , k B 

- 

10 2 , 10, 10 2 

W m _1 K _1 

coupling factor 

G 

- 

10 1 ® 

W m -3 K _1 

conductances 

C’e? O p 

- 

10 5 " 11 , 10 7-9 

W m” 2 K 1 

diffusivities 

Ote ) CUp j OL s 

CX.i = ki jCi 

10~ 2 , 10"®, 10" 4 

m 2 s 1 

conductance times 

Tee) Tcp 

T~c i = Cih/oi 

10 4 - 10“ 2 , 10 4 - 10 2 

ps 

diffusion times 

e? 7"dp 

~ T~di ~~ h j OLi 

1, 1000 

ps 

relaxation times 

Te, T p 

n = Ci/G 

1, 100 

ps 

Biot numbers 

B e , B p 

Bi =L — hai/ki 

10 -4 - 10 2 , 0.1 - 10 

- 

M numbers 

M e , M p 

Mi = Tdi / Ti = Gh/oi 

1, 10 

- 

R number 

R 

R = - Tp j Te == Cp / Ce 

1000 

- 


TABLE III: Quantities used in this work. Subscripts i=e, p, or s refer to the metal electrons, the metal phonons and the 
substrate phonons respectively. 


1. Fast diffusion approximation 

To obtain compact formulas, we consider here the two-temperature model in the fast diffusion (FD) approximation: 
the temperature is assumed to be uniform within the film, implying that both electrons and phonons diffusion are 
neglected. With the vector notation T = [T e , T p ] t , where T denotes the transpose, the equations are 


d t T = — MT, 


The decay is found to be exponential, with a characteristic decay time r given by 


Tr - \/Tr 2 - 4D ’ 

Tr(M) = - + — + — + 


( 8 ) 


Tr 

D = Det(M) = 


T e T p 

1 

Te^cp 


' cp 


' p ' ce 


'ce'cp 


To get a qualitative understanding of the r dependence on a e , we first consider the two limiting values for a e —> 0 
and cr e -A oo, which are denoted as r max and r m i n . Physically, T max represents the metal’s cooling time in the absence 
of any electron-phonon conductance, while r m , n is the cooling time in the limit of high electron-phonon conductance. 
On physical grounds, one can expect that upon increasing a e , r decreases monotonically from r max to T m ; n . To 
simplify even further, we use the fact that r e is often small compared to all others time involved in the problem, i.e. 
r e < Tp,T ce ,T C p. Then 


' cp 


1 + - + 0(r e 2 ) 


TnTd 


+ 0(r e ). 


(9) 

( 10 ) 


' cp 


Compared to Eq. [8j both expressions are numerically accurate and no further term in the expansion is necessary. As 
regards the ratio r, the above approximations yields 

TVnin 1 


(1 + T e /r p )(l + t cp /t p ) ’ 
1 


1 + T cp / 


for r e < r p . 


This shows that at low cr p values, r depend largely on a e . For instance, taking er p = 10' [44| yields r ~ 1/100, i.e. 
one obtains a 100-fold reduction in r when a e increases from 0 to oo. 






















For conductances a p smaller than 10 8 , r is actually rather well approximated by the simple expression 


t-t„ 


7m ax TVr 


1 + <T e /a p ' 


( 11 ) 


Deviations amount to only a few percents for a p < 10 s , and occur only at very large a p , where a p should be replaced 
with a function <?(<r p ) which increases slightly less than a p (see below). The resulting cooling time is shown in Fig.0 
which demonstrates a strong effect of the electron-phonon interfacial conductance <r e , especially for low values of cr p . 
A useful approximation for the following may be obtained, if we consider that r m ; n — X P +° P — T p and r max ~ r cp , 
which is a reasonable approximation given that r cp r p with the default parameters. This yields 


' cp 


1 + (Te/cr p ’ 


( 12 ) 


2. Full two-temperature model 

The FD model neglects diffusion of the energy carriers in the metal, which may induce additional resistance to the 
heat flow. Here, we go one step further as compared with the FD model, and account explicitely for the electron and 
phonon diffusion. As detailed in Appendix A, the long-time behavior of the temperature decay depends on whether 
the substrate temperature is constant or not. In the latter case, the decay is power-law, without any characteristic 
time scale. In the former case, the decay may be characterized by a single decay time r. This is the reason why we 
use this cold substrate approximation and determine analytically how the characteristic decay time r depends on cr p 
and cr e . We find that if a e ^> a p , the temperature decay is much faster than that expected for zero a e . 

Using h and r e as the length and time units, the relevant equations, boundary and initial conditions read as: 


d t T e = ^ dl x T e - (T e - T p ), 
Rd t T p = dl x T p + (T e - T p ), 


d x T e 

= 0 , 

d x T p 

= 0 , 

for 

x = 0. 

d x T e 

= B e T ei 

-d x T p 

— B p T p , 

for 

x = 1 

T e 

= T\ 

T P 

= 0, 

for 

t =0, 


There are thus six dimensionless parameters: B e , B p , M e , M Pl R and S/h. The general solution of the previous set of 
equations is not exponential. However, in order to characterize the long time decay, we consider Laplace transforms 
with respect to time and take its expansion close to zero frequency. As further explained in the appendices A and B, 
this gives an estimate of the longest relaxation time characterizing the long time cooling kinetics. Even if this relies 
on an approximation, we believe that it can reflect faithfully the influence of a e . 

For the two-temperature model considered here, the characteristic decay time r can be obtained analytically using a 
symbolic computation software, but the expression is much too large to be reported. As shown in Fig.[5]for the default 
parameters and a wide range of a p and cr e , large a e can result in a significant decrease of the decay time. Alternately, 
this trend may be captured analytically as detailed in the Appendix B. Finally, as visible in Fig. [5j the FD model 
provides a very reasonable approximation to the full two-temperature model, showing that even when “neglecting” 
diffusion, one can recover the qualitative behavior. This suggests that the effective conductance characterizing the 
different energy channels present at a metal-non metal interface, does not depend primarily on the diffusion of the 
energy carriers in the metal. 


C. Effective thermal boundary conductance 

From the solution of the two models analyzed above, we can predict the effect of the different electron-phonon 
processes on the effective thermal boundary conductance er e ff- This latter quantity is the apparent interfacial conduc¬ 
tance extracted from a TTR signal using a one-temperature model. Because the one-temperature model ignores the 
electronic degrees of freedom, the effective conductance depends on both the electron-phonon coupling constant and 
the conductance a e . Here, we discuss this dependence. The effective conductance is computed using 

tTeff = q/{T p - T 0 ), (13) 

cLT 

where q = ~h-^- is the interfacial heat flux density. Note that here we are assuming the phonon temperature to be 
uniform, in the spirit of the FD model as we concluded that phonon diffusion within the metal film does not change 
significantly the metal cooling kinetics. 
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FIG. 5: Characteristic decay time r as a function of the electron-phonon conductance <j e , for several phonon-phonon conductance 
<jp, as computed with the FD model and with the full two-temperature (dashed and solid lines respectively). 


Because in the FD model discussed previously, the long time decay of the metal lattice temperature is exponential, 
with a time constant r the effective conductance Eq (fT3l) is given by: er e ff = c p h/r. Using Eq. (THU) for the characteristic 
time, we predict the following expression of the conductance: 

l/<Teff = l/Gh + l/(cr e + cr p ). (14) 

This equation has a simple interpretation: the effect of the electron-phonon processes is equivalent to add in series 
the bulk resistance l/Gh with the total interfacial resistance l/(cr e + <7 p ). This latter resistance comprises two parallel 
interfacial resistances quantifying electron-phonon and phonon-phonon processes at the interface. At this point, it is 
important to realize that the expression of the bulk resistance l/Gh differs from the expression 1/y/Gk p derived by 
Majumdar and Reddy 0- The difference stems from the different boundary conditions considered for the electrons: 
Ref. 0 assumed an insulating boundary condition for the electrons, as they were not allowed to directly exchange 
energy with the substrate. Here, because we consider the possibility of a direct exchange at the interface, the boundary 
condition is not adiabatic. This difference of boundary condition changes the expression of the internal resistance 
governed by electron-phonon processes in the metal thin film. 

Finally, we discuss the traditional separation between equilibrium and non-equilibrium contributions to the effective 
conductance 0, [jOI- The expression of the conductance Eq. [HI has been derived under non-equilibrium conditions. 
If we had assumed equilibrium, we would have obtained the intuitive expression of the interfacial conductance <r e + c p 
which expresses the fact that the two energy channels -electrons and phonons- act in parallel to transmit heat to the 
substrate. Hence, we find that the non-equilibrium contribution consistent with an isothermal boundary condition is 

O-neq = Gh, (15) 

and is found to be larger for thick metal films h > 100 nm, and metals having a large electron-phonon coupling 
G > 10 17 W m -3 K _1 . For these systems, the non-equilibrium conductance is a neq > 10 10 W m -2 K _1 and the 
corresponding resistance is practically negligible, as compared with the phonon-phonon and electron-phonon interfacial 
resistances. On the other hand, for thin Au films with h < 10 nm the non-equilibrium conductance is (7 neq < 250 
MW m -2 K _1 and may become comparable with the interfacial conductances er p and a e . We conclude then that in 
metal-non metal superlattices with a small period, non-equilibrium effects can not be neglected, and should contribute 
to enhance the thermal resistance. 

We now discuss the effect of a finite a e on the value of the effective conductance Eq. (TTH) . All the expressions of 
the characteristic time derived so far clearly indicate that the electron-phonon conductance a e has a strong effect on 
the long time decay of the temperature of the metal, and may significantly accelerate the cooling kinetics. To further 
illustrate this point, we have considered the relative enhancement of the metal cooling induced by the electron-phonon 
conductance. This enhancement is simply defined by the ratio x = < T eff(o'e)/c r eff(o'e = 0) and as such may be considered 
as a simple estimate of the increase of the apparent thermal boundary conductance induced by a e . Figure [6] shows that 
the relative enhancement may reach one or two orders of magnitude, depending on the value of the phonon-phonon 
conductance cr p . The enhancement is most significant when the electron-phonon conductance a e becomes greater 
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than (jp, a situation that might be common, as discussed in Sec. [TT] Again, we conclude that the largest enhancements 
are observed when the metal/substrate interface transmits poorly the phonons, which corresponds to the case of 
acoustically mismatched solids. 



FIG. 6: Relative enhancement x °f the apparent thermal boundary conductance due to the electron-phonon conductance a e , 
as predicted by the fast diffusion approximation Eq. [14] with default parameters. 


IV. TWO REALISTIC CASES: GOLD/SILICON AND GOLD/SILICA INTERFACES 

A. Methods and parameters 


We have discussed above the effect of a e using a set of typical values for the different transport and thermodynamic 
quantities characterizing the metal and the dielectric. Here, we consider two realistic cases, which allows to take 
into account the temperature dependence of the transport coefficients. We focus on the gold/silicon and gold/silica 
systems, for which the values of the electron-phonon interfacial conductance are known [231 ] 



k 

C p 


-?Debye 

n 

CT e 

CTp 

Gold phonons 

18.0 

2.35 x 

To®" 

165 

97970 

- 

- 

Gold electrons 300.0 

1.97 x 

10 4 

- 

97970 

- 

- 

Silicon Si 

153.6 

1.68 x 

10 6 

625 

82891 

134 

72.6 

Silica Si02 

1.370 

1.01 x 

10 6 

403 

44075 

152 

141.5 


TABLE IV: Thermophysical parameters at 300 K: thermal conductivity k (Wm 1 K 1 ), heat capacity C p (Jm 3 K 1 ), Debye 
temperature Tbebye (K), molar density n (molm~ 3 ), electron-phonon cr e and phonon conductance <r p (MWm _ 2 K _1 ). 

The electron-phonon coupling factor was set to the value G = 2.5 x 10 10 MW m ” 3 K _1 predicted by DFT calculations 
when the electronic temperature is below 5000 K [28j. The values of the other thermophysical parameters at 300 K 
are given in Tab. IIVI The heat capacity of the gold electrons is calculated with C e = yT e where 7 = 65.64 J m -3 K ” 2 
is Sommerfeld’s constant of gold, while the vibrational heat capacity is obtained from standard expression [35}. 
As regards the thermal conductivity of the substrate, we used experimental values for silicon and silica extracted 
from Refs. [36] and [[37] respectively, as briefly described in the Appendix C. Finally, we take the electronic thermal 
conductivity and the phonon thermal conductivity to be constant in gold metal. 

The phonon-phonon thermal boundary conductances are assumed to be described by the diffuse mismatch model, 
as discussed in Sec. EH Because gold metal has a low Debye temperature @d = 165 K, we can assume that above 200 
K, the phonon-phonon conductance is constant, and takes the value 



where n is the number density, and the phonon transmission coefficient a has been calculated in eq. [ 6 ] Numerically, 
this yields cr p = 72.6 MW m -2 K _1 for the Au/Si interface, and er p = 141.5 MWm -2 K _1 for the Au/Si 02 interface. 
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The temperature dependence of the electron-phonon conductance for silica and silicon was extracted from Ref. 
using the following expressions 


CT Au/Si = 25.18 + 0.363 T, (17) 

CT Au/Si0 2 = 96.12 + 0.189 T, (18) 

where temperature is given in K, and conductances in MWm -2 K -1 . Equations [171 and [151 provide a reasonable 
description of the experimental data, as shown in Fig. [3] 

The equations of the two-temperature model are solved using a finite-difference scheme, with perfectly matched 
layers at the boundaries so as to simulate a quasi-infinite substrate, whose thickness is fixed at 50 fin i. Gradients are 
estimated using finite differences with three successive lattice points, while the condition h(T(x + ,t)—T(x~,t)) = kd x T 
is enforced at the interface using two points to evaluate the local gradient. The typical lattice step is 5 nm; the timestep 
is chosen so that the maximal temperature increment between two successive time steps is AT = max x \T(x,t n +i) — 
T(x,t n )\ = 0.01 K if t < 1 ps; AT = 0.5 K otherwise unless the temperature difference becomes |T — To| < 2 K, when 
a finer criterion AT = 0.05 K is retained. Here and in all the following To is the initial temperature of the system, 
prior to any heating. The initial value of the electronic temperature is T e = To + 2000 K, the metal penetration 
depth is S = 20 nm and the film thickness is fixed to h = 100 nm. 


B. Results 

We present now the results obtained by numerically solving the two-temperature model, using the physical parame¬ 
ters detailed above. As in Sec. cm we focus on the electronic temperature of the metal at the air/metal surface, since 
this is a simple measure of the metal reflectivity, which is probed in the time-resolved thermoreflectance experiments. 

Figure [7] a displays the temporal evolution of the surface electron temperature for the gold/silicon system. For 
each temperature, curves with the electron-phonon conductance switched off ( a e = 0) are also shown. The difference 
between the two cases is striking. Clearly, discarding the electron-phonon channel results in an electron cooling which 
is much slower. After 1000 ps, the electronic temperature is almost twice larger in the case of zero a e . Figure 0a 
thus illustrates the key point of this work: in the time interval relevant to thermoreflectance measurements 100— 1000 
ps, the electron-phonon interfacial coupling has a strong effect on the temperature decay. Therefore it should be 
taken into account, along with the phonon-phonon conductance. It should be acknowledged at this point, that for the 
nanosecond time regime we are considering here, electrons and phonons in the metal are practically at equilibrium as 
illustrated in Fig. [8] A slight difference between electron and phonon temperature is however visible when a e = 0, as 
a result of the adiabatic boundary condition obeyed by the electrons in this case. 

To quantify the effect of a e on the temperature decay, we introduce a time ti/ 2 defined as the time after which 
the electronic temperature has cooled down to half its initial value [45|. t\/ 2 is a simple estimate characterizing the 
cooling kinetics, that may be used in experiments. Figure [9] displays T\/ 2 as a function of the initial temperature To. 
The effect of the interfacial electron-phonon coupling is clearly seen: in presence of this channel, the electron cooling 
is typically a factor two to three faster than in its absence. In terms of apparent thermal boundary conductance, this 
corresponds to a change by an amount ranging from 200 to 300%. Again, we see that interfacial electron-phonon 
coupling should contribute significantly to the cooling process and to the apparent value of the thermal boundary 
conductance. Note also that Fig. [5] shows cooling times that are practically independent of the initial temperature 
To, in both cases where a e has been set to zero or not. This may certainly be explained by the use of a constant 
phonon-phonon conductance in the model. The use of a temperature dependent electron-phonon conductance a e does 
not change the dependence on To. 

We have also characterized the effect of cr e on the cooling kinetics of the gold/silica system. As illustrated in Fig. [3b, 
the temperature decay is weakly dependent on the interfacial electron-phonon conductance. This is confirmed by Fig. [9] 
where T\/ 2 is seen to be practically insensitive to the presence of the additional interfacial energy channel. This relative 
difference between the cases of gold/silicon and gold/silica may be explained by two factors. The first is the difference 
in relative magnitude of the two conductances: while a e is twice larger than a p for gold/silicon, they are nearly equal 
in the case of gold/silica. The other factor is the difference of thermal conductivities of the substrates, with silicon 
having a conductivity two orders of magnitude higher than silica. As a result, any energy flux coming from the metal’s 
electrons is conveyed away from the interface. In contrast, heat in the silica remains confined in the vicinity of the 
interface, slowing down any energy transfer coming from the metal’s electrons, and resulting in the relative insensivity 
to <7 e in the cooling kinetics of the gold/silica system. 

As a final illustration of the relevance of the electron-phonon boundary conductance, we consider for the case of a 
gold/silicon interface a conductance <r e varying in a wide range. The time T\/ 2 is shown in Fig. 1101 where the abscissa 
axis has been normalized by the ’experimental’ value given by Eq. 1171 The decrease of the cooling time ti/ 2 is very 
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FIG. 7: (Color online) Effect of the electron-phonon interfacial conductance <j e on the cooling of gold metal’s electrons, as 
found with the two-temperature model for gold film on silicon substrate (a) and on silica substrate (b). Solid lines: simulations 
with (Je given by Eas. 1171 and 1181 : Dashed lines: same simulations with u e = 0. 


sharp for low <r e , suggesting that even a moderate electron/substrate coupling may lead to an enhanced thermal 
conductance. For larger values of cr e , the decrease is somewhat less pronounced, but still leading to a threefold faster 
cooling rate for the highest value studied. Again, this roughly corresponds to a factor three in the apparent thermal 
boundary resistance. Note that the corresponding value of a e ~ 1300 MW m -2 K _1 may be realistic to describe metals 
other than gold, as we discussed in Sec. |TT] 


V. CONCLUSION 

In summary, we used a combination of analytical work and numerical simulations of the two-temperature model to 
probe the influence of a direct interfacial electron-phonon energy transfer in the thermal transport between a metal, 
or a semi-metal and a dielectric substrate. We have shown that Sergeev’s model predicts for most metal-non metal 
systems considered here an electron-phonon conductance cr e which is larger or comparable with the phonon-phonon 
conductance a p . Perhaps surprisingly, even semi-metals such as Bismuth, with poor electronic transport properties 
compared to those of noble metals, may involve a significant <r e , but which remains smaller than er p . This suggests 
that electrons of materials which are poor electronic conductors may participate in interfacial heat transport. 

To assess the effect of a finite electron-phonon interfacial conductance, we investigated the cooling of a metal thin film 
instantaneously heated by a strong laser pulse. This situation is of direct relevance to time-resolved thermoreflectance 
technique, which is used to measure the thermal boundary resistance between the metal and the substrate. The 
effective interfacial conductance consists in a non-equilibrium contribution Gh and an equilibrium conductance cr p +cr e . 
The non-equilibrium conductance is found to be important for metals having a relatively low value of G, and for 
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FIG. 8: (Color online) Electron and phonon temperature profiles across the gold/silicon interface, 1000 ps after initial heating. 
Solid lines correspond to simulations with <j e given by Eq. 1171 while dashed lines are simulations without interfacial electron- 
phonon coupling <r e = 0. The temperature profiles of the substrate are also shown. 
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FIG. 9: (Color online) Cooling time t \/2 (see text) for the gold/silicon and gold/silica systems, as a function of the temperature 
To prior to heating. Open symbols are simulations where a e given by eas. 1 171 and 1181 while crosses and stars correspond to the 
case a e = 0. 


thin films a situation relevant to metal-dielectric superlattices or core-shell nanoparticles. The effect of a e is found to 
enhance the apparent thermal boundary conductance as compared with <t p , up to two orders of magnitude for acoustic 
mismatched interfaces, characterized by low values of a p on the order of 10 MW m -2 K _1 . This may explain the large 
discrepancies reported between the apparent conductances measured for solids having a large Debye temperature ratio 
and the classical DMM model @. Numerical simulations of the two-temperature model on two realistic cases confirm 
our analytical findings: while a e has only a mild effect on the cooling of gold/silica system, it significantly speeds 
up the cooling of the gold/silicon system. Those distinct behaviors originate in the difference in the phonon-phonon 
conductance cr p and the thermal conductivities of the substrate. 

On the theoretical side, one final remark is in order. We emphasized the strong limitations of Sergeev’s model to 
predict the electron-phonon conductance a e . In particular, Sergeev’s model disregards the nature of the substrate, 
in disagreement with the experimental results [23l |. and seems to overestimate the interfacial electron-phonon con¬ 
ductance, in the few cases when comparison with experiments is possible. Hence, there is a clear need to develop 
theoretical models for the electron-phonon conductance that would account for the nature of the non-metal substrate. 
Given the small number of a e values that have been obtained experimentally so far, measurements for a variety of 
metal/substrate systems would be invaluable in this endeavour [271 ]. 
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FIG. 10: Effect of the electron-phonon conductance a e on the cooling time t\i 2 (see text) for the gold/silicon system. Here cf e 
denotes the value extracted from experiments [2Sj and given by Eg. 1171 The temperature is To = 300 K. The line is a guide to 
the eye. 
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Appendix A: Characterizing the long-time temperature decay. 


This appendix focuses on the long-time behavior of the temperature decay. We find that when the substrate heating 
is negligible, which we call the “cold substrate” approximation, the long-time decay may be characterized by a single 
decay time r. On the other hand, if the substrate heating is significant the long-time decay is power-law, with no 
characteristic time. This distinction is essential when one wants to extract the thermal boundary conductance from 
the experimental curves. For simplicity, we use the simple one-temperature model to make this point, but it applies 
to the two-temperature model as well. 

Cold substrate. Here we assume that the temperature increase of the substrate is negligible, that is T s = 0 
at all time. Taking h as the unit length and the phonon conductance time T cp = c p h/a p as the unit time, the 
adimensionalized equations for the one-temperature model are 

B p d t T p = d 2 x T p , 

d x T p =0 for x = 0, 

B P T P + d x T p =0 for x = 1, 

T p = T i {x) for t =0. 


There are only two parameters: the Biot number B p = ha p /k p and the S/h ratio. We now introduce T p (x,s ), the 
Laplace transform of T p with respect to time, which satisfies the equation 


whose solution is 


-d 2 xx f p + B p fp = B p T\x), 


T p {x,s) = Cie 7 * + C 2 e“ 7X + Tpart > 7 = \/B^s, 

T P art = — f sinh [y(u — x)] T l {u)du. 

7 Jo 


The constants C\ and Ci are obtained from the boundary conditions at x = 0 and x = 1. From now on, we focus on 
the temperature at the air-metal interface T p (s) = T p (0,s), and consider the long-time behavior by taking a small-s 
expansion. The result is T p (s) ~ ao + ais + 0(s 2 ). Note that for a time dependent function /(£), whose Laplace 
transform has the expansion /(s) = ao + a\S + 0(s 2 ), we can always characterize the long time decay by an effective 
time r e ff = — ai/ao. If fit) ~ e -t / T is a single exponential, one recovers r e g- = r. If on the other hand, /(£) is a sum 
of exponential terms: 


then the procedure yields 


T p {t) =Y J c ne~ t/T \ 

n 



(x) = 


J2 C n 


(19) 


7~eff — 


( 20 ) 
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that is an average that takes into account the amplitude of each mode, and gives more weight to the slow modes. 

Following this procedure for the metal phonon temperature, we obtain the characteristic time r for arbitrary 
penetration depth S 


3(B P S - 1) {B p (2S 2 - 1) - 2) - 2ei (B p (3B P <5 3 - 3(B p + 1)5 2 + B p + 3) + 3) 
6(e*(B p ($-l)-l)-B p <y+l) 

—2e h ^ (k 2 (h 3 — 3 hS 2 + 3d 3 ) + 3 n(h — S)(h + 6) + 3/i) — 3 (5k — 1) ( h 2 n + 2 h — 2 5 2 k) 

6 a p K ( e h / s (h,K — 5n + 1) + Sk — l) 


( 21 ) 

( 22 ) 


Equation [21] applies if length and time units are h and r cp respectively. Equation [22] is given in terms of dimensional 
parameters, and for simplicity of notations, we have introduced the inverse length k = a p /k p . This formula simplifies 
in the two limiting cases 



T 

Bl 

(23) 

<5 —>■ 0, 

Tcp 

= 1 + 

3{B p + 1) 

5 —> oo, 

r 

1 + -B p (5B P + 8) 

(24) 


Tcp 

12(B p + 2) 


The Biot number B p indicates the relative importance between the heat flux within the bulk of the film and the flux 
at the interface. If B p <C 1, the temperature is quasi-uniform within the film and the flux is limited by the interface 
resistance, then r = r cp , whatever the initial profile. On the other hand, if B p 1, diffusion within the film is the 
limiting factor and r ~ h 2 /a p — r cp is the typical time of phonon diffusion over the film thickness h , with a prefactor 
that depends on the initial profile. Both limiting cases are visible in Fig. 1111 which shows the characteristic time r 
as a function of a p . Interestingly, with the default parameters, a typical value <r p = 10 8 yields B p = 1, meaning that 
none of the limits applies. And indeed, the difference is significant between r cp = 1000 ps and r =1167, 1256 and 
1361 ps for 5 = 0, 20, oo respectively. 



FIG. 11: Characteristic decay time r in the one-temperature model as a function of the conductance ct p , with default parameters, 
and penetration depth <5 = 0, 20 nm and oo, from bottom to top. 


Heated substrate. Retracing the steps of the analysis above, one finds for the small-s expansion of the surface 
temperature 


(25) 


When taking the inverse Laplace transform, the terms with half-integer exponents yield power-laws, which are domi¬ 
nant over other terms [4(| . We thus expect for T p ( t ) the approximation 


T P (t) 



0 - 1/2 

2v / tt 


3 

t 2 + 


3 03/2 s 
4i/7T 


+ o(H) 


(26) 


The key point is that the decay is power-law, since at the longest time, T p (t) — yfwt. This is a general feature 

of diffusion on a half-line [38| . Furthermore, the coefficient a_i /2 = hc p /y/c s k s depends on the total energy injected 
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and on the substrate properties, but not on the interfacial conductance a p [47J. In contrast with the cold substrate 


approximation treated above, it is not possible to define a single time scale that could characterize the temperature 
decay. Thus, those calculations with the one-temperature model show that two cases should be distinguished. If 
substrate heating is negligible, the decay may be characterized by an effective time r, and presumably the decay time 
r can be easily extracted, from which, if all other parameters are known, o p can be deduced. If on the other hand, the 
cold substrate approximation does not apply, the temperature decay is power-law, and at the longest time, completely 
independent of cr p . It might be possible, using fits of the experimental curve, to estimate the coefficients a p / 2 , 03/21 
and subsequent, and deduce a p , but the process seems much more delicate. This conclusion applies a fortiori to the 
two-temperature model. 


Appendix B: Two-temperature model in the cold substrate approximation. 


Here, we apply to the full two-temperature model 
the dimensionless parameters 


the analysis detailed in Appendix A. To do so, let’s introduce 


a 2 — M e + M p , 
h2 = (M 2 + RM 2 ) 2 
4(M e + M p )3 ’ 
2 M e M p (1 + R ) 
Me + Mp 


(27) 

(28) 

(29) 


The full expression for the characteristic time r is very unwieldy, in part because r involves power of e a , an exponential 
term can not be easily approximated since a is neither small nor large: with the default parameters, a = 3.31 is of 
order unity. 

To get a qualitative understanding of the r dependence on <r e , we first consider the two limiting values for a e —> 0 
and <r e —> 00, which are denoted as r max and r m i n . Even after taking the limits, the expressions for T max and r m ; n 
are still untractable. We therefore resort to two bold approximations. First, the limit 8/h —»■ 0 is taken, even if the 
value 0.2 is not so small. Second, a large-a expansion is made, even though e a ~ 27 is not so large. The resulting 
expressions are quite simple 


Tmax 

T e 

Tmin 

T e 


1 + ir. 


w + a-i+0[- 


c 2 ; 6 b + c 
3 " + 3a 


36 + 2c 2 
3 a 2 


+ 0 



and yields for the ratio r = r m ; n 


B p M e 

3a 2 



1 



(30) 

(31) 


(32) 


Our main conclusion follows from those approximate expressions. Large values of cr e can induce, compared to the 
a e = 0 case, a very pronounced decrease in r. For instance, tr p = 10 7 yields r = 0.012, that is a 100-fold reduction. 


Op 

Tmax 

App. 

FD 

Tmin 

App. 

FD 

10' 

10288 

10314 

10100 

90 

104 

99 

5 10 7 

2213 

2234 

2020 

88 

- 

95 

10 8 

1208 

1224 

1010 

86 

- 

91 

10 9 

322 

315 

101 

78 

- 

50 


TABLE V: Maximum and minimum decay times r max and r m i n . App. and FD correspond to Eqs. J30J and © for the former, 
and Eqs. m and <EB for the latter. The conductance o p is expressed in W m 2 K 1 and the times in ps. 


Table IVl compares exact numerical values for r m ; n , r max , and their approximation given by Eqs. (l3Ul) - (l3Tl) . which 
are surprisingly accurate, given the assumption on S and a. Note that the approximate Eq. m predicts that r m ; n 
is independent of <r p . However, when taking the five-term expansion of r m ; n , one recovers the slight decrease with er p 
which is observed for the exact r m j n . Finally, one can obtain a good approximation of r with the simple formula 


r — T m in _ f 

'T — T • 1-p _^_ ’ 

'max ' mm r T 


( 33 ) 
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FIG. 12: The function g(a p ) as defined from Eg. 1331 (points!. The line indicates the function g(a p ) = a p . 


where g{cr p ) is a function plotted in Fig. [T2] For low cr p < 5 x 10 7 , g(a p ) ~ a p while it becomes smaller at higher a p . 

Two side remarks are in order. First, the formulas given above for the coefficients a, b and c are valid for any 
dimensionless parameters but may be simplified in many cases. For instance, taking the default parameters of 
Tab. IIII1 yields M e = 1, M p = 10, R = 1000, which suggests taking the limiting case R 3> 1, and M e -C M p , and 
leads to the simple expressions 


a 1 2 3 4 5 6 7 8 9 10 11 12 = M p , b 2 = c 2 = M e R. 

Second, there is a simple limiting case for r max as given by Eq. [30] In addition to R 1, let’s assume that r max r e 
and B p -C 1. With default parameters, the former approximation is well justified and the latter applies only for low 
<jp values (B p = 0.1 for fr p = 10' for instance). Then the first term inside the brackets of Eq. (13U1) dominates over the 
others, and one finds r max = r cp , i.e. one recovers the characteristic time for the one-temperature model in the limit 
B p —> 0, as given in Appendix A. 


Appendix C: thermal conductivity for silicon and silica. 

We have used the following fits to describe the temperature dependence of the thermal conductivity: 

k S i = 203913 T -1 - 26 , (34) 

fc Si o 2 = 0.624lnT-2.19, (35) 

where the temperature T is expressed in K and the conductivities in W m -1 K -1 . The experimental values and fitting 
curves are displayed in Fig. [T3l 
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